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High-mass X-ray binary systems are powered by the stellar wind of their donor stars. The X-ray 
state of Cygnus X-l is correlated with the properties of the wind which defines the environment of 
mass accretion. Chandra-HETGS observations close to orbital phase 0 allow for an analysis of the 
photoionzed stellar wind at high resolution, but because of the strong variability due to soft X-ray 
absorption dips, simultaneous multi-satellite observations are required to track and understand 
the continuum, too. Besides an earlier joint Chandra and RXTE observation, we present first 
results from a recent campaign which represents the best broad-band spectrum of CygX-1 ever 
achieved: On 2008 April 18/19 we observed this source with XMM-Newton , Chandra , Suzaku, 
RXTE, INTEGRAL, Swift, and AGILE in X- and y-rays, as well as with VLA in the radio. 
After superior conjunction of the black hole, we detect soft X-ray absorption dips likely due to 
clumps in the focused wind covering > 95 % of the X-ray source, with column densities likely to 
be of several 10 23 cm~ 2 , which also affect photon energies above 20 keV via Compton scattering. 
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Introduction 

One of the currently still puzzling issues in our understanding of microquasars is the exis- 
tence of the different states (low/hard, high/soft, very high/steep power-law, . . . ), which distinguish 
themselves by the soft X-ray luminosity, the X-ray spectrum, the timing properties, and also the 
nature of the jet and its radio spectrum. One assumes that these states with their in many ways 
different properties correspond to different configurations of the accretion flow which cause dif- 
ferent channels or efficiencies of energy dissipation. For a high-mass X-ray binary system, the 
global framework of mass accretion onto the black hole is defined by the powering wind of the 
donor star. One possibility is that changes in the wind properties (like density, velocity, ionization 
parameter, . . . ) might finally trigger state transitions. 

This contribution focuses on the wind of the well-known persistently bright X-ray source 
CygnusX-1. Its companion star HDE226868, a ~18M 0 09.7 type supergiant [11], fills 97% 
of its Roche lobe [6]. The stellar wind is therefore strongly focused due to tidal interactions [5,7]. 
The detailed structure of this wind defines the environment of mass accretion onto the black hole on 
the largest scale. At a luminosity of To ,5-250keV = 3 x 10 37 erg s 1 (which is obtained already in the 
low/hard state [9]), the X-ray source produces a considerable feedback on the wind by photoion- 
ization [2], which can create a complex wind structure within the binary system. The parameters 
inferred from the Ha emission line are (anti-)correlated with the X-ray flux: in the low/hard state, 
the wind is probably denser and therefore less strongly ionized, which allows for a more efficient 
wind acceleration and thus a lower accretion rate onto the black hole [8], 

We have obtained multi-satellite observations of Cyg X- 1 in the hard state close to superior 
conjunction of the black hole, i.e., at phase d orb « 0 in the 5.6 day binary orbit [3, 8], which arc 
ideal to investigate the wind: at the low inclination i ~ 35° [11], no eclipses of CygX-1 by its 
companion star can be observed. Nevertheless, absorption dips, which arc likely related to the 
focused stream covering the X-ray source, arc often detected close to d orb = 0 [1], 

1. The Chandra and RXTE observation of 2003 April 
1.1 Advantages of a multi-satellite observation 

In the first paid of this contribution, we show some of the results from a simultaneous obser- 
vation of Cyg X- 1 performed by Chandra and RXTE in 2003 April shortly before 0 or b = 0. (More 
details arc given by Hanke et. al [9, 10].) The light curves shown in Fig. 1 clearly show the ad- 
vantages of a multi-satellite observation-. Due to RATE’S low Earth orbit, there are data gaps of 
~45 minutes length in each RXTE-orbit (lasting ~90 minutes) from passages through the South 
Atlantic Anomaly (SAA). The high variability which is present on timescales of hours and which 
is related to the stellar wind could hardly be investigated with the RXTE data alone. The (continu- 
ous) soft X-ray light curve taken by Chandra is shaped by absorption dips. Only the joint spectral 
analysis of the non-dip continuum allows to combine the strength of the different instruments. For 
example, the PCA spectrum clearly requires an Fe Ka emission line at 6.4 keV, but the resolution 
is not good enough to decide whether it is narrow or relativistically broadened. In the simultaneous 
Chandra spectrum, a narrow component is detected (see Fig. 2), but its flux is not sufficient to ex- 
plain the PCA feature, from which we infer the presence of a broad iron line as well. (More details 
on iron lines in Cyg X-l arc given by M. A. Nowak in his contribution to this volume [16].) 
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Figure 1: Light curves of CygX-1 during the simultaneous Chandra/RXTE observation of 2003 
April. Top: RXTE-HEXTE (20-250keV). Middle: RXTE-PCA (4-20 keV). Bottom: Chandra-HETGS 
(0.5-12keV 1 st order events). The dashed lines define the non-dip times. 


1.2 The non-dip spectrum 

The non-dip spectrum can be obtained by filtering for high count rates, see the dashed horizon- 
tal line in Fig. 1. The results of this subsection are thoroughly described by Hanke et al. (2008) [9]. 

The Chandra -HETGS spectrum of the bright source CygX-1 shows neutral absorption sig- 
natures of the interstellar medium (and probably also source environment). We detect structured 
absorption edges, namely the Fe L 2 and L 3 edges at 17.2 and 17.5 A (0.721 and 0.708 keV) - 
which J. C. Lee describes in detail in her contribution in this volume -, the O K-edge at 22.8 A 
(0.544 keV) with the Ka resonance absorption line, and the Ne K-edge at 14.3 A (0.867 keV) with 
the K/3 resonance absorption line. These complex absorption edges can be modelled with tbnew 1 , 
an improved version of tbvarabs, which contains those fine structures. 

The non-dip spectrum is dominated by absorption lines from the stellar wind of highly ion- 
ized ions, namely H- and He-like O, Ne, Na, Mg, Al, Si, S, Ar, Ca, and Fe, as well as many Fe 
L-shell absorption lines, see Fig. 2. We fit these features with Gaussian absorption lines and infer 
equivalent widths and Doppler shifts with respect to the rest wavelengths of the identified transi- 
tions. In a second approach, we have defined a fit-function (in ISIS [12, 15]) for the whole series 
of absorption lines of an ion. Theoretical Voigt profiles were used which take into account the 
oscillator strengths of the strongest transitions. The ion’s column density is a direct fit-parameter 
of the model, avoiding an explicit curve of growth analysis. The velocity shifts inferred from both 
methods are rather low , mostly |v| < 200km s _1 . The latter result can be understood within a sim- 
ple model for the spherical wind, e.g., a generalized CAK-model [4, 14]: 


v(r) = v(r) • - with v(r) = vq + (voo-vq)- 1 


R * 


See http : / /pulsar . sternwarte . uni-erlangen . de/wi lms/ resear ch/tbabs. 


(1.1) 
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Figure 2: Cfiandra-HETGS CygX-1 non-dip spectrum with identification and fits of absorption lines of 
H- and He-like ions (Lyman and Helium series; black labels) and Fe L-shell ions (blue labels). Note the 
narrow and broad Fe K a line (red) at 1.94 A (=6.4 keV, see second a- (energy-) axis on top of each panel). 
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Figure 3: Wind velocity, projected against the 
black hole. Absorption lines with low velocity 
v ra d can be observed from the dark purple re- 
gion close to the black Thales circle of v' ra( j = 0, 
v ra d > 0 can only be observed from the red re- 
gion inside this circle, and v ra d < 0 only from 
the blue region outside of it. The brighter the 
red or blue, the higher is |v ra d(r)|. The dashed 
gray lines are contours of constant v ra d in steps 
of 200 km s“ 1 . The highlighted area with the p orb 
scale shows the region in the orbital plane with 
| cos $ | < sin /, which equivalently corresponds 
to all possible lines of sight (for i = 35°). 


For the velocity observed in X-ray absorption lines, one has to consider the projection against 
the black hole , v ra d(r) = cos a(f) ■ v(r) , where a(r ) is the angle between wind velocity v(r) and 
the line of sight towards the black hole. Figure 3 shows v ra d (r) as a colorscale plot. On the cir- 
cle passing through the center of the star and the black hole, a = 90° by Thales’ theorem, and 
thus v ra d = 0 - independent of the assumed velocity law, as long as it is spherically symmetric. 
The circle separates (possible) absorbing regions of the wind seen at a redshift (inside) from those 
seen at a blueshift (outside). The location of absorbers seen at a low velocity is thus well con- 
strained. It is important to realize that this result does not depend on the model for the wind 
velocity, provided that v(r) points radially away from the star. 

Due to the system’s inclination, not every line of sight in the orbital plane can actually be 
observed: the angle # between the inclined line of sight l = (sin/, 0, cos/) and the binary axis 
a = (cos(27 r0 ol -b), sin(27T0 orb ), 0) is given by cos # = l o a = sin/ • cos(27T0 orb ). The range of 
possible $ values is obviously limited by a low inclination. (For a face -on geometry with / = 0°, 
there is clearly no orbital variation, i.e. # = 90°, while every 0 < if < 180° is possible for edge-on 
systems with i = 90°. ) The region with | cos #| < sin / is highlighted in Fig. 3 between the crossing 
solid gray lines. For a spherically symmetrical wind, this region in the orbital plane is equivalent to 
all lines of sight experienced by an inclined observer. The green labels show the position of these 
lines of sight for given orbital phases 0 orb , where d orb and 1 — </) ol - b are identical. 

As a result, the observation of low Doppler velocities in absorption lines in front of the X-ray 
source constrains the location of the absorber - rather independently of the wind model. This 
location is in good agreement with XSTAR simulations of the photoionization zone. The number 
density of Fl-atoms, n H , is obtained from the velocity profile (Eq. 1.1) and the continuity equation, 
A/* = /./ «?h ■ /i h (r) -4 nr 2 ■ v(r). Flere, p = 1.4 is the mean molecular weight per Fl-atom and 
»7 h is the mass of an Fl-atom. Inserting the numbers = 2 1 00 km s 1 and R+ = 11 Rq [11] gives: 

x -2 

«h(t> = 2.2x 10 10 cm~ 3 -d(r/Rf) , where d{x) = w (1.2) 

vo/voc + (1 -vo/voc) • (1 - 1/xY 

Between the stellar surface and the black hole, x = r/R, varies between 1 and 2.4 and the average 
wind density is /in ~ 10 11 cm -3 for vq/voo = v t herm/v<>° = 0.01 and a reasonable j8 « 1. 
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Figure 4: The Mg xi triplet of 
resonance (r) absorption and in- 
tercombination (i) as well as for- 
bidden (f) emission lines is de- 
tected with two pairs of i- and 
f-lines in the non-dip Chan- 
dra-HETGS spectra of Cyg X-l. 
The red-shifted pair stems from 
an emitting plasma component 
of higher density and can be 
identified with the focused wind. 


While these previous findings from the non-dip spectrum could be explained by a spheri- 
cal wind only, there are also signatures of a focused wind : Figure 4 shows the Mg XI triplet 
of ls 2 ( 1 So) — > I s2p (' 1 P i ) “resonance” (r), ls 2 ( 1 So) <— ls2p( 3 Pi.2) “intercombination” (i) and 
Is 2 ('So) <— 1s2s( 3 Si) “forbidden” (f) transitions. The i- and f-emission lines show up as two 
components each: one pair is almost at rest - in accord with the r absorption line while the other 
pair is seen at a strong redshift (500-900 km s 1 ). The ratio R = F(f)/F(\) of the line fluxes is usu- 
ally used for density diagnostics [17]. However, in the presence of strong UV-radiation fields (such 
as in the vicinity of HDE 226868), the R ratio is systematically decreased due to photodepletion of 
the 1s2s( 3 Si) level [13]. Therefore, an absolute density analysis would overestimate the density. 
One can, however, still infer the relative statement that the red-shifted pair of i- and f-lines (with 
the lower A’-ratio) stems from a much denser emitting plasma component. Both its velocity and its 
density indicate that this component can be identified with the focused wind. 

1.3 The dip-spectrum 

During the dips (Fig. 5a), the soft X-ray spectrum is considerably absorbed. Defining the 
soft X-ray color as the ratio of the 0.5-1. 5 keV and the 1.5-3 keV Chandra count rate, and the 
hard color as the ratio of 1.5-3 keV and 3-10keV rate, a color-color diagram (Fig. 5b) shows that 
neutral photoabsoiption alone cannot explain the data. The nose-like shape of the data points can 
be explained by a partial covering model for the photoelectric absoiption: 

(1 - /c) • + f c . e -Ah,2-<Tbf(£) (1.3) 

We use a constant = 5.4 x 10 21 cm~ 2 and a power-law continuum with photon index F = 1.73 
as found with RXTE-FCA for the non-dip spectrum [9]. In this model, the strongest dips can be 
described by an additional absorber with > 4 x 10 23 cnF 2 covering a fraction / c between 95 % 
and 98 % of the X-ray source. 

The high-resolution dip spectra (see, e.g., Fig. 6) show strong absorption lines of lower ionized 
ions, e.g., resonance K a absorption lines of Si VI-XI between 6.8 and 7.1 A, which are not detected 
in the non-dip spectrum at all. This result is consistent with an origin of the dips in colder matter 
in the focused wind which covers the X-ray source. (More details will be presented in [10].) 
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Figure 5: (a): Chandra-HETGS (0.5-10keV) light curve (see also Fig. 1). (b): Color-color diagram from 

(0.5-1.5keV):(1.5-3 keV) and (1.5-3 keV):(3-10keV) count rate ratios in 25 s intervals, explained by a grid 
of partial covering models (Eq. 1.3), assuming a power-law continuum (r = 1.73). The gray labels indicate 
the covering fraction f c in % and the column density (Vh , 2 of the additional absorber in 10 22 cm 2 . 



Figure 6: Chandra -HETGS CygX-1 spectrum of the first dip ( 0 orh < 0.96, see Fig. 5a). While the non-dip 
spectrum (Fig. 2) showed silicon resonance K a absorption lines of only Si XIV and XIII (Si Fya at 6.18 A 
and Si Hea at 6.65 A), the dip spectrum also presents the lower ionization stages Si VI-XI (green labels). 
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Table 1: Multi-satellite observation of CygnusX-1 on 

2008 April 18/19. 

Satellite 

Start 

Stop 

Notes 

XMM 

INTEGRAL 

2008-04-18, 12:07:27 
2008-04-18, 12:35:32 

2008-04-19, 04:25:46 
2008-04-19, 23:00:45 

C jobs, of the Cyg-region} 

RXTE | 
Suzaku 

2008-04-18, 13:21:36 
2008-04-19, 17:37:20 
2008-04-18, 16:22:33 

2008-04-19, 04:22:39 
2008-04-19, 20:04:47 
2008-04-19, 10:08:08 

j- with SAA-interruptions 
with SAA-interruptions 

Chandra | 

2008-04-18, 18:27:19 
2008-04-19, 14:57:27 

2008-04-19, 02:48:58 
2008-04-19, 20:17:26 



Note that, according to [8], 0 or b(2OO8-O4-18, 15:53) = 0 and 0 or b (2008-04- 19, 20:00) = 0.20 . 


2. The multi-satellite campaign of 2008 April 

For a holistic study of CygX-1 and the focused wind of HDE 226868, further multi-satellite 
observations had to be performed. The strong variability due to dipping requires to observe Cyg X- 1 
strictly simultaneously, if one wants to benefit from several different instruments in order to study 
the wind and the absorption dips. We originally proposed for a joint XMM+Chandra 2 observation 
- with XMM’s EPIC-pn camera for a high sensitivity for the spectral continuum shape through 
the broad iron line region, Chandra's HETGS for high-resolution wind spectroscopy, and XMM's 
RGS for a high-resolution spectrum at energies below the O-edge. RXTE was available to provide 
the broadband continuum spectrum up to 250 keV. We were also able to obtain simultaneous INTE- 
GRAL time - extending the spectral coverage up to ~2MeV as well as Suzaku - allowing for 
an independent confirmation of the spectral continuum, see Table 1. We were furthermore awarded 
a few ks of simultaneous Swift TOO time. AGILE did not detect the source in the 30 MeV-50 GeV 
band. Cyg X- 1 has also been monitored in the radio by the VLA during this observation. 

As expected, the observation was again dominated by dipping events. Most of them arc re- 
stricted to the soft X-ray band ( Chandra , XMM), but a few arc also detected above 4 keV by the 
PCA on RXTE or even above 20 keV by IBIS on INTEGRAL. A color-color diagram reveals the 
same nose-like shape as was found for the dips in 2003 (Fig. 5b), indicating the same possible 
origin of the dips in partial covering of the X-ray source [16, 10]. 

Due to the fast decrease of the photoelectric absorption cross section o^f above the ionization 
threshold. X-ray photons with energies well above the iron edge at 7. 1 keV can ha idly be absorbed 
by the wind. They can, however, still be scattered out of the line of sight. Figure 7 shows such an 
X-ray scattering trough lasting more than 8 hours: Between April 18, ~ 1 6:45 and April 19, ~ 1 : 1 5, 
the 20-40 keV count rate measured with INTEGRAL -IBIS, as well as the 12-60 keV count rate 
measured with Suzaku-PIN, decrease until April 18, 20:15 and rise again afterwards. A continuous 
model joining two linear functions (a descending one followed by an ascending one) can describe 
both light curves simultaneously (with a scaling factor of 0.23 between the PIN and IBIS count 
rates) remarkably well (^ r 2 ed = 2). If the 30 % reduction between the highest and the lowest count 
rate of this fit (corresponding to non-dip and the strongest dip) is due to Compton scattering, the 
optical depth T = 0.36 corresponds to /V H = 6x 10 23 cm~ 2 . 

2 Note that scheduling an simultaneous observation with orbital phase constraints is not easy for these two satellites... 
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Figure 7: Light curves obtained with INTEGRAL -IBIS in the 20-40 keV band (top black data, left y axis) 
and Suzaku -PIN in the 1 2-60 keV band (lower blue data, right y-axis) with a time resolution of 200 and 96 s, 
revealing an X-ray scattering trough shortly after superior conjunction, which lasts for 8.5 hours. During 
this time, the light curves can be described by two continuously jointed linear functions (red lines). 
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Figure 8: The 0.5keV-2MeV broadband spectrum of CygX-1, observed simultaneously with XMM, 
Suzaku, RXTE and INTEGRAL in 2008 April. The data below 300 keV can be described consistently by 
a broken power-law with exponential cutoff at high energies (Lfoid ~ 170keV). ISIS’ (model independent) 
flux-corrected spectra (relying on the response matrices) seem to indicate a non-thermal tail in the SPI data. 

The broadband spectrum between 0.5 keV and 2 MeV, assembled from various instruments 
(Fig. 8), reveals that CygX-1 was again observed in the hard state. Below 300 keV, it can be 
described as a simple elbow-shaped broken power-law with high energy cutoff. RXTE-PCA’s 
photon indices below and above Ebreak = 9.8keV, namely f| = 1 .7 and 44 = 1.4, are amongst the 
hardest ones found in the long term monitoring of Cyg X-l from 1999 to 2004 [16, 18]. Preliminary 
INTEGRAL -SPI data seem to indicate a non-thermal tail in the spectrum above 400 keV, but further 
analysis through joint broadband fits with Comptonization models is still ongoing. 
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